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Abstract 

Social dynamic opinion models have been widely studied to understand how interactions among 
individuals cause opinions to evolve. Most opinion models that utilize spin interaction models 
usually produce a consensus steady state in which only one opinion exists. Because in reality 
different opinions usually coexist, we focus on non-consensus opinion models in which above a 
certain threshold two opinions coexist in a stable relationship. We revisit and extend the non- 
consensus opinion (NCO) model introduced by Shao et aZ.[l|. The NCO model in random networks 
displays a second order phase transition that belongs to regular mean field percolation and is 
characterized by the appearance (above a certain threshold) of a large spanning cluster of the 
minority opinion. We generalize the NCO model by adding a weight factor W to individual's 
own opinion when determining its future opinion (NCOVK model). We find that as W increases 
the minority opinion holders tend to form stable clusters with a smaller initial minority fraction 
compared to the NCO model. We also revisit another non-consensus opinion model based on 
the NCO model, the inflexible contrarian opinion (ICO) model |2 ; ], which introduces inflexible 
contrarians to model a competition between two opinions in the steady state. Inflexible contrarians 
are individuals that never change their own opinion but may influence opinions of others. To place 
the inflexible contarians in the ICO model we use two different strategies, random placement and 
one in which high-degree nodes are targeted. In both strategies, the inflexible contrarians effectively 
decrease the size of the largest cluster of the rival opinion but the effect is more pronounced under 
the targeted method. All of the above models have previously been explored in terms of a single 
network. However human communities are rarely isolated, instead are usually interconnected. 
Because opinions propagate not only within single networks but also between networks, and because 
the rules of opinion formation within a network may differ from those between networks, we 
study here the opinion dynamics in coupled networks. Each network represents a social group or 
community and the interdependent links joining individuals from different networks may be social 
ties that are unusually strong, e.g., married couples. We apply the non-consensus opinion (NCO) 
rule on each individual network and the global majority rule on interdependent pairs such that two 
interdependent agents with different opinions will, due to the influence of mass media, follow the 
majority opinion of the entire population. The opinion interactions within each network and the 
interdependent links across networks interlace periodically until a steady state is reached. We find 
that the interdependent links effectively force the system from a second order phase transition, 
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which is characteristic of the NCO model on a single network, to a hybrid phase transition, i.e., 
a mix of second-order and abrupt jump-like transitions that ultimately becomes, as we increase 
the percentage of interdependent agents, a pure abrupt transition. We conclude that for the 
NCO model on coupled networks, interactions through interdependent links could push the non- 
consensus opinion type model to a consensus opinion type model, which mimics the reality that 
increased mass communication causes people to hold opinions that are increasingly similar. We 
also find that the effect of interdependent links is more pronounced in interdependent scale free 
networks than in interdependent Erdos Renyi networks. 

PACS numbers: 
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I. INTRODUCTION 



Statistical physics methods have been successfully applied to understand the cooperative 
behavior of complex interactions between microscopic entities at a macroscopic level. In 
recent decades many research fields, such as biology, ecology, economics, and sociology, 
have used concepts and tools from statistical mechanics to better understand the collective 
behavior of different systems either in individual scientific fields or in some combination of 
interdisciplinary fields. Recently the application of statistical physics to social phenomena, 
and opinion dynamics in particular, has attracted the attention of an increasing number 
of physicists. Statistical physics can be used to explore an important question in opinion 
dynamics: how can interactions between individuals create order in a situation that is 
initially disordered? Order in this social science context means agreement, and disorder 
means disagreement. The transition from a disordered state to a macroscopic ordered state 
is a familiar territory in traditional statistical physics, and tools such as Ising spin models 
are often used to explore this kind of transition. Another significant aspect present in 
social dynamics is the topology of the substrate in which a process evolves. This topology 
describes the relationships between individuals by identifying, e.g., friendship pairs and 
interaction frequencies. Researchers have mapped the topology of social connections onto 
complex networks in which the nodes represent agents and the links represent the interactions 
between agents [3 15 1 . Various versions of opinion models based on spin models have been 



proposed and studied, such as the Sznajd model [Iff, the voter model [17|, |l8[, the majority 



rule model 
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20j, and the social impact model [21 




Almost all spin-like opinion models mentioned above are based on short-range interactions 
that reach an ordered steady state, with a consensus opinion that can be described as a 
consensus opinion model. However, in real life different opinions are mostly present and 
coexist. In a presidential election in a country with two political parties in which each party 
has its own candidate, for example, a majority opinion and a minority opinion coexist. The 
opinions among the voters differ, with one fraction of the voters supporting one candidate 
and the rest supporting the other, and rarely will the two opinions reach consensus. This 
reality has motivated scientists to explore opinion models that are more realistic, ones in 
which two opinions can stably coexist. Shao et al. |l| proposed a nonconsensus opinion 
(NCO) model that achieves a steady state with two opinions coexisting. Unlike the majority 
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rule model and the voter model in which the dynamic of an agent's opinion is not influenced 
by the agent's own current opinion but only by its neighbors, the NCO model assumes that 
during the opinion formation process an agent's opinion is influenced by both its own current 
opinion and the opinions of friends, modeled as nearest-neighbors in a network. This NCO 
model begins with a disordered state with a fraction / of cr + opinion and a fraction 1 — / 
of cr_ opinion distributed randomly on the nodes of a network. Through interactions the 
two opinions compete and reach a non-consensus stable state with clusters of o~ + and ct_ 
opinions. In the NCO model, at each time step each node adopts the majority opinion of 
its "neighborhood" , which consists of the node's nearest neighbors and itself. When there is 
a tie, the node does not change its opinion. The NCO model takes each node's own current 
opinion into consideration, and this is a critical condition for reaching a nonconsensus steady 
state. Beginning with a random initial condition, this novel nontrivial stable state in which 
both majority and minority opinions coexist is achieved after a relatively short sequence of 
time steps in the dynamic process. The NCO model has a smooth phase transition with 
the control parameter /. Below a critical threshold f c , only the majority opinion exists. 
Above f c , minorities can form large spanning clusters across the total population of size 
N. Using simulations, Shao et al. |l| suggested that the smooth phase transition in the 
NCO model in random networks is of the same universality class as regular mean field 
(MF) percolation. But simulations of the NCO model in Euclidean lattices suggest that the 
orocess might belong to the universality class of invasion percolation with trapping (TIP) 



23(| . Apparently this is the first time, to the best of our knowledge, that a social dynamic 



model has been mapped to percolation, an important tool in statistical physics. However, 
the nature of this percolation on 2D lattice is still under debate 23j, |24( . Exact solutions of 
the NCO model in one dimension and in a Cayley tree have been developed by Ben-Avraham 



25|. 

Here we present simulations suggesting that the behavior of the NCO model, in which two 
opinions coexist, disappears when the average network degree increases. When the average 
degree of a network is high, the agent's own opinion becomes less effective and the NCO 
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model converges to the majority voter model. This was argued analytically by Roca 
and claimed also by Sattari et al [24]. In the present paper, we also generalize the NCO 
model and create a nonconsensus opinion model by adding a weight (NCOiy model) to an 
agent's own opinion. The weight W > 1 represents the strength of an individual's own 
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opinion. Note that in the NCO model W is assumed to be 1 like the weight of its neighbors 
opinion. We find that the NCOW model inherits all the features of the NCO model, except 
that the critical threshold f c of the NCOPU model with W > 1 decreases when W increases. 
This means that strengthening one's own opinion helps smaller minority opinion groups to 
survive. 

The NCO model reaches a steady state in which the two opinions coexist. This is only 
partially realistic. In real life, two opinions do not simply coexist — they continue to com- 
pete. Real-world examples include the decades-long competition between the Windows and 
Macintosh operating systems and between Republicans and Democrats in US presidential 
politics. All the participants in these competitions have the same goal: winning. In order 
to increase their prospects of winning, they need as many supporters (or customers) as pos- 
sible. Thus, it is interesting to study how two opinions continue to compete after they have 
reached a steady state. In order to consider both aspects, the nonconsensus steady state 
and the competition, Li et al. [2] proposed an inflexible contrarian opinion (ICO) model 
in which a fraction <p of inflexible contrarians are introduced into the final steady state of 
the NCO model and two different competition strategies are then applied. The concepts 



28| in his work 



of inflexible agents and contrarian agents were introduced by Galam 
on opinion models. In the ICO model, an inflexible contrarian is an agent that holds an 
opinion contrary to that held by the majority of its surrounding group and its opinion is 
not influenced by its surrounding group — it never changes. Inflexible contrarians have one 
goal: to change the opinion of the current supporters in the rival group. We see this strat- 
egy when, for example, companies send a free product to potential customers in order to 
convince them to adopt the product and influence their friends to do the same. We study 
the ICO model in order to determine, for example, whether these free products actually do 
help to win the competition, how many free products are needed to be sent, and who are the 
best candidates to receive the free product. Reference Q presents two strategies for intro- 
ducing inflexible contrarians into the steady state of the a+ opinion groups: (i) the random 
strategy and (ii) the targeted (high degree) strategy. Using these strategies, we find that 
the relative size of the largest cluster in state a + undergoes a second-order phase transition 
at a critical fraction of inflexible contrarians <p c below which the two opinions can coexist 
and above which only a_ exists. Thus the ICO also belongs to the type of nonconsensus 
opinion models. The results also indicate that the largest cluster in state a + undergoes a 
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second order phase transition that can be mapped into MF percolation similar to the NCO 
model. 

All opinion models described above have been studied on a single network. However, 
in real social opinion dynamics, individuals belonging to different social communities can 
as well communicate. In a traditional agrarian village, for example, two separate working 
relationship networks often form. Men work in the fields with other men and women work 
in their homes with other women. Marriages between men and women in this setting create 
interdependencies between the two working relationship networks. As far as we know, there 
has been no model study of how this kind of strong social connection between two such 
different groups influences the exchange of opinions. In studying the opinion dynamics across 
different groups we utilize a concept that has recent ly g ained wide attention: the resilience 



of interdependent networks to cascading failures |29l-l36l|. Connecting two networks together 
with interdependent links allows individuals to exchange opinions between networks. In 
our model, two nodes from different networks that are connected by interdependent links 
represent a pair of nodes that have strong social relations. In interdependent networks we 
usually distinguish between the connectivity links between agents within each network or 
community and interdependent links between agents from different networks. 

To study the effect of interdependent links on opinion dynamics, we propose a non- 
consensus opinion (NCO) model on coupled networks in which we assume different opinion 
formation rules for internal connectivity and interdependent links. We assume that during 
the dynamic process of opinion formation the agents that are connected with interdependent 
links will have the same opinion, this being the case because their social relationship is 
strong. In our model, the NCO rules are applied in each individual network. For the 
coupled pairs the following rule is applied: if two interdependent nodes have the same 
opinion, they will keep this opinion, but if they have different opinions, they will follow the 
majority opinion of the interdependent network system (global majority rule). Many other 
possible rules could be tested for the interdependent pairs, but we adopt here, for simplicity, 
the majority rule. When an opinion is shared by two interdependent individuals, such as 
a married couple, because their social relationship is strong and close, they will tend to 
maintain their opinion against outside influence. If their opinions differ initially, they tend 
to eventually resolve their differences and share the same opinion. In the process of resolving 
their differences, however, they can be significantly influenced by outside forces, e.g., mass 
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media, and thus we assume that they often end up sharing the majority opinion. When we 
increase the number of interdependent links between the coupled networks, the transition 
changes from a pure second order phase transition to a hybrid phase transition and finally 
to a seemingly abrupt transition. The hybrid transition contains both a second order and 
an abrupt transition. The model type of the NCO model on coupled networks also changes 
as the number of interdependent links increases, and thus the system goes from being a kind 
of nonconsensus opinion model to being a kind of consensus opinion model. This suggests 
that strong interactions between different social groups is pushing our world in the direction 
of becoming more uniform in their opinions. 

The paper is organized as follows, in Sec. [TT] we revisit some important concepts on the 
topology of opinion clusters and percolation. We then present the results and discussions 
on NCO and NCOPF model in Sec. ED on ICO model in Sec. |TV] and on NCO model on 
coupled networks in Sec. [V] Finally, we present our summary in Sec. [V]] 



II. TOPOLOGY OF OPINION CLUSTERS AND PERCOLATION 



In recent decades, many researchers have studied how network topology affects the pro- 



cesses that evolve in them. Examp 



diseases, and percolation |4H6|. 
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es of such processes are the spreading of rumors, opinions, 
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43]. Classical percolation processes deal with the ran- 



dom failure of nodes (or links) and present a geometrical second order phase transition with 
a control parameter p that represents the fraction of nodes (or links) remaining after a 
random failure of a fraction 1 — p of nodes (or links). 

There exists a critical probability p c above which a "giant component" (GC) appears. 
The number of nodes in the GC, Si, is called the order parameter of the phase transition. 
Below criticality there is no GC and only finite clusters exist. For p < p c the size distribution 
of the clusters is n s ~ s~ T with a cutoff that diverges when approaching p c . At criticality, 
in the thermodynamic limit, the size of the second largest component S% diverges at p c 
as S 2 ~ \p — Pc|~ 7 just as the susceptibility with the distance to the critical temperature. 
For large networks (N — > 00), p c = 1/(k — 1), where k is the branching factor given by 
K = (k 2 ) I (k), where (k) and (k 2 ) are the first and second moments of the degree distribution 
P(k) of the network respectively ll|. We perform all our simulations on both Erdos-Renyi 
(ER) networks and scale-free (SF) networks [h]]. ER networks are characterized by 
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a Poisson degree distribution, P(k) = e~^ k > (k) k /k\. In SF networks the degree distribution 
is given by a power law, P(k) ~ k~ x , for A; min < k < k max , where k m i Q is the lowest degree 
of the network and & max is the highest degree of the network. For random SF networks 
^max ~ A rl /( A ~ 1 ) is the degree cutoff [lit, where N is the system size and A is the broadness 
of the distribution. 

We begin by examining percolation in ER networks. At criticality, percolation in ER 
networks is equivalent to percolation on a Cayley tree or percolation at the upper critical 
dimension d c = 6 where all the exponents have mean field (MF) values with r = 5/2 and 
7 = 1. Note that in the ER case, the mass of the incipient infinite cluster S± scales as iV 2//3 at 
criticality. We can understand this result by using the framework of percolation theory for 
the upper critical dimension d c — 6. Since Si ~ R df and N ~ R d (where d is the dimension 
of the initial lattice, df the fractal dimension, and R the spatial diameter of the cluster), it 
follows that Si ~ iV d // dc and since d c = 6 and d f = 4 we obtain Si ~ iV 2/3 [id]. 

For SF networks, the GC at criticality is Si ~ iV 2/3 for A > 4, and Si ~ j\r(^2)/(A-i) for 
44( | . For SF networks, with A < 3, (k 2 ) — > oo when N — > oo because L M -» oo 



3 < A < 4 



and thus p c = 0, making these networks extremely robust against random failures 

However if we decrease fc max by targeting and removing the highest degree nodes (hubs), 
p c is finite 45] and we recover a second order phase transition with MF exponents as for 
ER networks. We will show below that this is also true for our model here. A similar MF 

n 

behavior in SF networks with A < 3 was found also by Valdez et al. [46|] for the percolation 
of susceptible clusters during the spread of an epidemic. In our simulations we always choose 
kmm = 2 for SF networks in order to ensure that they are almost fully connected [11]. 



III. THE NCO MODEL 

In the NCO model [l| on a single network with N nodes, opinion a + and cr_ are initially 
randomly assigned to each node with a fraction of / and 1 — / respectively. The basic 
assumption of the NCO model is that the opinion of an agent is influenced by both its 
own opinion and the opinions of its nearest neighbors (the agent's friends). The opinion 
formation rule states that at each time step, each node adopts the majority opinion, which 
includes both the opinions of its neighbors and itself. If there is a tie, the node's opinion 
will remain unchanged. Using this rule, each node is tested at each simulation step to see 



9 



whether its opinion has changed. All these updates are performed simultaneously and in 
parallel until no more changes occur and a steady state is reached. 

Figure [T] demonstrates the dynamic behavior of the NCO model on a small network with 
nine nodes. At time t — 0, five nodes are randomly assigned opinion a + (empty circle), and 
the remaining four, opinion cr_ (solid circle). After checking the status of each node, we find 
that only node A belongs to a local minority with opinion cr+, so at the end of this time step, 
node A changes its opinion to cr_. At time t=l only node B belongs to a local minority, so 
at the end of this time step, the opinion of node B will be updated to er_. At time t=3, 
every node has the same opinion as its local majority, where the final nonconsensus steady 
state is reached. 

A. Simulation Results 

In the steady state s± = S\/N is the normalized size of the largest opinion a + cluster, 
S2 = S2/N is the normalized size of the second largest opinion a + cluster, and F is the 
normalized fraction of opinion cr + nodes. Figure [2] shows plots of s±, S2, and F as a function 
of the initial fraction / of the opinion a + nodes for both ER and SF networks. We find that 
due to the symmetrical status of both opinions, F is a monotonically increasing function of / 
with symmetry around / = 0.5. Figure [2] also shows the emergence of a second order phase 
transition. Note that there is a critical threshold f c , which is characterized by the sharp 
peak of S2- Below f c , s± approaches zero, where only the majority opinion can form steady 
clusters, and above f c , s± increases as / increases and a state with stable coexistence of both 
majority and minority opinion clusters appears. Although for both ER and SF networks 
f c < 0.5 as expected, ER networks have smaller values of f c than SF networks for the same 
average degree (k). For example for SF networks with fc min = 2 and A = 2.5 where (k) ~ 5.5, 
f c ~ 0.45, while for the same average degree (k) = 5.5, for ER networks, f c ~ 0.4. These 
differences in f c indicate that the minorities in SF networks need more initial supports to 
form final steady state clusters, than minorities in ER networks. This might be understood 
due to the high degree nodes (hubs) of a SF network. In the NCO model a hub, because 
of its large number of connections, is strongly influenced by its neighbors which are with 
high probability of the majority opinion, and puts the minority opinion at a disadvantage. 
Reference [1( presents also studies of the NCO model in a two-dimensional Euclidean lattice 
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and of the NCO model on real-world networks. 

We next present numerical simulations indicating that the phase transition observed in 
the NCO model is in the same universality class as regular MF percolation. Percolation in 
random networks (e.g., ER and SF networks with A > 4 for random failures or all A for 



targeted attacks) 13j, |42j, |43) is obtained by MF theory, which predicts that at criticality the 



cluster size distribution is n s ~ s~ T with r = 2.5 and S± ~ N e , where 9 = df/d c with df = 4 
and d c = 6 represent the fractal and the upper critical dimension of percolation respectively 
and thus 9 = 2/3 (See Sec. [TTJ) . Figure [3](a) shows the finite cluster size distribution n s of the 
a + opinion cluster as a function of s at criticality (/ = f c ). Figure MJo) show Si at criticality 
f c as a function of N for ER and SF networks with (k) = 4 and A = 2.5, respectively. Note 
that in both networks r ~ 2.5 and 9 w 2/3. These two exponents strongly indicate that the 
NCO model in random networks behaves like a second order phase transition that belongs 
to the same universality class as regular MF percolation. 

In our above results we focus on networks that have a relatively low average degree (k) . 
We test the model for networks with higher average degrees. In SF networks we increase 
(k) by increasing the value of fc m i n . Figure H] shows si and S2 as a function of / for different 
values of (k) for ER networks and SF networks. As the values of (k) increase, i.e., as the 
network becomes increasingly condensed and the number of interactions between agents 
increases, a sharper change of Si at a critical threshold is observed. This may suggest (but 
can not be proved by simulations) the existence of a critical value (k) = k c that is strongly 
affected by the topology of the network. Below k c , as (k) increases, f c shifts to the right, as 
can be seen from the shift of the peak of S2- Above f c two opinions can continue to coexist 
and remain stable. Above k c the smooth second order phase transition is replaced by a 
sharp jump of s± at approximately / = 0.5 that is accompanied by the disappearance of the 
peaks of S2- Note also that as the values of (k) increase, the region in which two opinions 
coexist becomes increasingly smaller and approach zero for very large values of (k) possibly 
above k c . In terms of the NCO model, as the number of connections between individuals 
increase, the opinion of each individual becomes less important and each individual becomes 
increasingly susceptible to the influence of the majority opinion across the entire system. 
Thus the majority opinion can easily overwhelm the minority opinion, causing the critical 
behavior of the NCO model, the second-order phase transition, to disappear at large (k) and 
the NCO model to converge to the majority voter model yielding a possible global consensus 
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throughout the system. Note that analytical arguments for an abrupt transition of the NCO 
model at large (k) are given in Ref. 261 ] . 

How can one help the minority opinions to survive? As we have seen, as the number 
of friends of an agent increases, the importance of the agent's own opinion decreases. In 
this way the majority opinion gradually eliminates the minority opinion. If we generalize 
the NCO model by adding a weight value W to each agent's own opinion, as W of an 
agent increases, the influence of the opinion of the agent's neighbors decreases. We call this 
generalization of the NCO model the NCOW model. As in the NCO model, in the ~NCOW 
model we change an agent's opinion if he is in a local minority but we also weight the agent's 
own opinion W times more than its nearest neighbors. The NCO model is actually a special 
case of the NCOiy model in which W — 1. Figure shows plots of Si and S2 as a function of 
/ for both W — 1 and W = 4. Note that as W increases, the second-order phase transition 
becomes flatter and the peak of the S2 shifts to the left, which indicates a smaller critical 
threshold f c . The smaller value of f c for larger values of W means that from the minority 
point of view, it needs fewer initial supporters to form and maintain stable finite clusters. 
When weight is added to the agents own opinion (indicating stubbornness) they become less 
susceptible to outside influence. Thus in the NCOiy model the majority is aided when the 
agents make more friends, but the minority in turn is aided when the agents treat their own 
opinion as more important than their friends' opinions. 



IV. THE ICO MODEL 



The initial configuration of the inflexible contrarian opinion (ICO) model corresponds 
to the final steady state of the NCO model in which two opinions a + and <r_ coexist. At 
t = a fraction of inflexible contrarians of opinion cr_ are introduced into clusters of 
cr+ by replacing nodes of <7+. The inflexible contrarians are agents that hold a strong and 
unchangeable <r_ opinion, that theoretically could influence the a+ opinion of their neighbors 
as the system evolves with NCO dynamics. Because the opinion held by the inflexible 
contrarians is unchanging, they function as a quenched noise in the network. The system 
evolves according to NCO dynamics until a new steady state is reached. In this steady state 
the agents form clusters of two different opinions above a new threshold f c = f c {4>). Because 
the contrarians hold the cr_ opinion, the size of the a + clusters decreases as <fi increases. 
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Figure [6] demonstrates the dynamic of the ICO model. We use two different strategies to 
introduce a fraction of inflexible contrarians. In strategy I we chose the fraction of 
nodes with a + opinion at random. In strategy II the inflexible contrarians are chosen from 
the agents with cr + opinion in decreasing order of their connectivity. Strategy II is thus a 
targeted strategy. 

A. Simulation Results 

We present our simulation results for ER networks with (k) = 4 and N = 10 5 . For 
simulation results of ICO model for SF networks see Ref. [2(|. Figure [7] shows plots of s± and 
S2 as a function of / for different values of for strategies I and II, respectively. Note that 
the ICO model inherits some of the properties of the NCO model. This is the case because 
there is a smooth phase transition with a critical threshold f c , where f c is characterized by 
the sharp peak of S2- However, for the ICO model, f c is also a function of 0. Thus we denote 
the new f c in the ICO model by / c (0). We find that, as increases, the critical value f c (4>) 
increases, which means that the largest cluster composed of cr+ agents becomes less robust 
due to the increase in the number of inflexible contrarians of opinion cr_. Note also that for 
/ > f c {4>), si decreases as increases. Thus, we conclude that inflexible contrarians with 
opinion cr_ have two effects: (i) they increase the value of / c (0) and thus the a + opinion 
needs more initial support in order to survive, and (ii) they decrease the size of the largest 
a + opinion cluster at / > f c {4>)- Note also that in the ICO model when is large the largest 
<j + cluster is fully destroyed and the second-order phase transition is lost. This is probably 
due to the fact that when is large, minority groups do not have high degree nodes and 
thus their average connectivity becomes smaller than 1 and, as a consequence, will no longer 
be able to form stable clusters [7|. As expected (see Fig. [7J strategy II is more efficient in 
destroying the largest minority component. This is plausible because, when selecting the 
initial fraction of inflexible contrarians using a targeted strategy, almost all the inflexible 
contrarians will be in the largest initial a + cluster since this cluster includes most of the 
high degree nodes. Figure |8]test this hypothesis and shows at the final stage of the NCO the 
ratio F(k), which is the number of nodes within the GC of a + opinion with degree k divided 
by the total number of nodes of opinion a + with degree k in the entire network system, for 
different values of /. We find that for large values of k, F(k) — > 1. These results support 
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our previous hypothesis that almost all the high degree nodes belong to the largest cluster, 
and this explains why strategy II is more efficient than strategy I. 

We next test whether the ICO model undergoes a phase transition as a function of and 
what it its type. Figures M^&) andlHKc) show plots of S\ as a function of for different values 
of / for strategy I and strategy II, respectively. Figures MJo) (top) and [^d) (top) show plots 
of S2 as a function of for different values of / for strategy I and strategy II, respectively. 
We can see that in both strategies s 2 has a peak at = <j) c {f), which is a characteristic of 
a second-order phase transition. Figures Mjo) (bottom) and |9^d) (bottom) further support 
the presence of a second order phase transition by showing plots of the derivative of Si 
with respect to for different values of /. Note that there is an abrupt change with in 
Asi/A0 at the same position of the peak of s 2 , suggesting that the transition is of second 
order. We next show that the second order phase transition has the same exponents as MF 
percolation. Figure [Till plots the finite cluster size distribution of a + agents, n s as a function 
of s at / = / c (0), from where we obtain r = 5/2. From s 2 we also compute the exponent 
7 and obtain 7 « 1 (not shown). These two exponents indicate that the ICO model on 
random graphs belongs to the same universality class as MF percolation. 

V. THE NCO ON COUPLED NETWORKS MODEL 

Figure [TTJ demonstrates the dynamics of the NCO model on coupled networks. In coupled 
networks, the two networks represent two groups of people. The links within each network 
denote the relationships between nodes. For simplicity, we assume that the two networks 
have the same number of nodes N and the same degree distribution. We also assign these two 
networks the same initial opinion condition, i.e., in both networks there is initially a fraction 
/ of nodes holding the a + opinion, and a fraction 1 — / holding the <r_ opinion. To represent 
the strong social coupling between the two groups, we randomly choose a fraction q of the 
nodes from both networks to form qN pairs of one-to-one interdependent pairs regardless of 
their original opinions. At time t=0, in both networks we apply the same opinion formation 
rule, the NCO model, to decide whether an agent will change its opinion regardless of 
the interdependent links. This means that at this stage opinions propagate in each single 
network independently — as though the other network does not exist. All opinion updates 
are made simultaneously and in parallel. At t=l, if two nodes with an interdependent link 
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have the same opinion they keep that opinion. If they do not, they follow the majority 
opinion of the coupled networks (global majority rule). All interdependent agents update 
their opinions simultaneously at the end of this time step. We repeat these two steps until 
the system reaches a steady state. 

A. Simulation Results 

We perform simulations of the NCO on coupled networks where both of the interde- 
pendent networks are either ER networks with (k) = 4 or SF networks with /c min = 2 
and A = 2.5. For an initial fraction / of opinion a + and a fraction q of interdependent 
links, the NCO on coupled networks is simulated on 10 4 network realizations to explore how 
interdependent links affect opinion dynamics. 

1. NCO on Coupled ER Networks 

We first investigate si as a function of / for different values of q. Figure fl2Ta) shows that 
when q = 0, which corresponds to the NCO model on a single network, the system undergoes 
a second order phase transition with a critical threshold f c When q > 0, there are two 
regions < q < 0.5 and q > 0.5. For the region < q < 0.5, as in the NCO model on a 
single network, the second order phase transition still exists, but the critical value / = f c (q) 
is increasing with q. The value of f c (q) can be determined by the location of the peak of S2, 
which is shown in Fig. fl2T b). where we plot S2 as a function of / for different values of q. The 
inset of Fig. |T2Ta) shows a plot of f c (q) as a function of q. We find that the peak of s 2 shifts 
to the right for q < 0.5 as q increases, which means that f c (q) increases as q increases. This 
suggests that if we add more interdependent links between the two networks, the minority 
opinion will need a larger initial fraction in order to exist. In the region < q < 0.5 we also 
find that, unlike the NCO model on a single ER network, there is an abrupt change of si at 
/ = 0.5, indicating that in addition to the smooth second order phase transition at f c (q), 
there may also be a discontinuous transition at / = 0.5. Our results suggest that when 
< q < 0.5 the system may undergo a hybrid phase transition [471 ] . which is a mixture of 
both an abrupt and a second order phase transition. We also find that as q increases the 
discontinuity around / = 0.5 becomes more pronounced. Although the system possesses a 
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seemingly discontinuous phase transition for < q < 0.5 , the model itself is still a non- 
consensus model, i.e., when / is above the critical value f c {q) the two opinions coexist in a 
steady state. When q > 0.5 the smooth second order phase transition of s± disappears and 
is replaced by an abrupt transition at / = 0.5. When q > 0.5 the peak of S2 disappears, 
supporting the loss of the second order phase transition [see Fig. IT21f b)]. and the system 
undergoes a pure abrupt transition. This suggests that when interactions between networks 
are sufficiently strong the hybrid phase transition is replaced by a pure abrupt transition. 
For all values of q, the region where two opinions can coexist decreases as q increases, and 
the NCO coupled networks model moves at large q toward the consensus type opinion model. 
To further support our finding of the existence of a discontinuous transition when q > 0, 
in Fig. [T21 c) and its inset we plot respectively the jump of Si, ASi = 5*1(0.51) — Si (0.49), 
and ASi/N, around / = 0.5 as a function of the system size N for different values of q. 
The linear relationship between A Si and iV supports our assumption of the existence of a 
discontinuous transition around / = 0.5 for all values of q > 0. Note that, as the value of 
q increases, the value of ASi increases, which means that as we increase the value of q the 
abrupt transition becomes more pronounced. 

To further support our conclusions, we investigate the number of iterations (NOI), which 
is the number of time steps needed to reach the steady state, as a function of / for different 
values of q. Figure [T2T d) shows a plot of the NOI as a function of / for different values of q. 
As described in Ref . |48] , in a pure first order phase transition due to cascading failures the 
location of the peak of the NOI determines the critical threshold of the transition, which 
is the case for q > 0.5 in our model. Figure fT2T d) shows that there is only one peak for 
the NOI curves for q > 0.5 at / = 0.5, which is the position of the critical threshold of the 
abrupt transition. In the hybrid phase transition for q < 0.5, the relation between the peak 
of the NOI and the critical threshold is unclear because there are two critical thresholds, 
one for the discontinuous transition at / = 0.5 and the other for the second order phase 
transition at f c (q)- Figure ITST d) shows that when q < 0.5 the NOI exhibits two symmetric 
peaks. The inset of Fig. H2T d) shows a plot of the location of the left peaks of the NOI as 
a function of q. Comparing the insets in Fig. fT2T a) and Fig. IT2f d). we find that the curve 
of the peak locations of the NOI is always above the f c {q) curve, which suggests that for a 
hybrid phase transition the peak of the NOI is located between the critical thresholds of the 
second order phase transition and the abrupt transition. 
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Figure [13] shows a log-log plot of the NOI at / = 0.5 as a function of the system size 
iV for different values of q, and the inset of Fig. [13] shows the same in a log-linear plot. 
The accuracy of the simulations is such that we cannot distinguish the relationship between 
exponential and logarithmic. However, the increase of NOI with system size indicates that 
there is a real jump at approximately / = 0.5 rather than a finite size effect. This supports 
our previous conjecture that for all values of q > 0, the NCO on coupled networks exhibits 
an abrupt transition at / = 0.5. 

We next present results indicating that, when q < 0.5 and when / is close to f c {q), our 
model is in the same universality class as regular MF percolation, even though a discontinuity 
appears at larger /. For regular percolation on random grap hs at criticality, the cluster sizes 



follow a power law distribution, n s ~ s T with r = 2.5 |13l . |42| . |43| . Figure IT41 shows a plot 
of function of s for finite cr+ clusters at criticality, f c (q)- We see that for q < 0.5, 

r m 2.5, and for q > 0.5, the power law no longer holds. The exponent values we obtain 
strongly indicate that, for small value of q, the NCO model on coupled ER networks close 
to f c is in the same universality class as mean field percolation in random networks. The 
power law for the cluster size distribution at q > 0.5 disappears, so we conclude that the 
NCO coupled networks model changes the phase transition type as q increases from q < 0.5 
to q > 0.5. 



2. NCO on Coupled SF Networks 

Empirical studies show that many real-world social networks are not ER. They instead 
exhibit a SF degree distribution [lo| in which P(k) ~ k~ x and A characterize the broadness 
of the distribution. A feature of SF is the existence of hubs, i.e., very high degree nodes. 
These large hubs make the opinion dynamic processes in SF networks much more efficient 



than in ER networks 



49 



54]. 



Because of its large number of connections, a hub in the NCO model tends to follow 
the opinion of the majority and effectively influence the opinions of its neighbors. In a SF 
network the hubs help the majority dominate the minority, and thus the NCO model on 



a single SF network has a larger f c and exhi 



ER networks with the same average degree jl]. This is also the case in interdependent SF 
networks. Figures IToTa). IToTb). and [ToTc) depict s\, S2, and NOI as a functions of / for 



)its a much sharper jump around f c than in 
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different values of q, respectively. The results of the NCO model on coupled SF networks 
are similar to those of coupled ER networks, except that the region of the hybrid phase 
transition is much smaller in coupled SF networks. This is confirmed by the fact that the 
peak of S2 drops much faster for small q values and that the single peak of NOI shows up 
at smaller q values for coupled SF networks in contrast to the case of coupled ER networks. 
This indicates that the pure abrupt phase transition occurs at smaller q values in coupled SF 
networks compared to coupled ER networks, which suggests that in coupled SF networks a 
smaller number of interdependent agents are needed to achieve a consensus state compared 
to coupled ER networks. Figure [TH] shows a plot of n s as a function of s for finite a + clusters 
at criticality. Note that in SF networks when q < 0.1 the n s decays as a power law with 
t = 2.5, and when q > 0.1 the power law decay of n s no longer holds. This suggests that 
only for small values of q our NCO model on coupled SF networks is in the same universality 
class as regular MF percolation. Comparing Fig. [16] with Fig. [T2J we find that the power 
law decay disappears at smaller q values in coupled SF networks compared to coupled ER 
networks. This supports our hypothesis that interdependent links push the entire system 
to an abrupt phase transition more effectively in coupled SF networks than in coupled ER 
networks. 

In both coupled ER and SF networks, our non-consensus opinion second order phase 
transition model is transformed into a consensus opinion type abrupt transition model when 
the number of interdependent links is increased. This suggests that increasing the inter- 
actions between different groups in our world will push humanity to become increasingly 
homogeneous, i.e., interdependent pairs in the NCO coupled networks model helps the ma- 
jority opinions supporters to eliminate the minority opinion, making uniformity (consensus) 
a possible final result. 

VI. SUMMARY 

In this paper we revisit and extended the non-consensus opinion (NCO) model, introduced 
by Shao et al. [l||. We introduce the NCOiy model in which each node's opinion is given 
a weight W to represent the nodes' resistance to opinion changes. We find that in both 
the NCO and the NCOW models the size of the largest minority cluster with a + opinion 
undergoes a second order MF percolation transition in which the control parameter is /. The 
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NCOW model is more robust than the NCO model because the weighted nodes reinforce 
the largest <t+ minority cluster, and shift the critical value f c to values lower than those 
found in the NCO model. We also show that when the average network degree (k) in 
the NCO model is increased, the second order phase transition is replaced by an abrupt 
transition, making the NCO model converge to a consensus type opinion model. We also 
review another non-consensus opinion model, the ICO model |2j, which introduces into the 
system, using both random and targeted strategies, a fraction <fi of inflexible contrarians 
(which act as quenched noise). As <fi increases, both random and targeted strategies reduce 
the size of the largest a + cluster and, above a critical threshold = <f) c , the largest a + 
cluster disappears and the second order phase transition is also lost. The targeted strategy 
is more efficient in eliminating the largest a + cluster or decreasing its size. This is due to 
the fact that the contarians are introduced (targeted) mainly into the largest cluster, which 
contains most of the high degree nodes. Thus a smaller C value is needed to eliminate 
the largest cluster of minority in the targeted strategy compared to the random strategy. 
We also study an opinion model in which two interdependent networks are coupled by a 
fraction q of interdependent links. The internal dynamics within each network obey the 
NCO rules, but the cross-network interdependent nodes, when their opinions differ, obey 
the global majority rule. These interdependent links force the system from a second order 
phase transition, characteristic of the NCO model on a single network, to a hybrid phase 
transition, i.e., a mix of a second order transition and an abrupt transition. As the fraction of 
interdependent links increases, the system evolves to a pure abrupt phase transition. Above a 
certain value of q, which is strongly dependent on network topology, the interdependent link 
interactions push the non-consensus opinion model to a consensus opinion model. Because 
scale free networks have large hubs, the effect of interdependent links is more pronounced 
in interdependent scale free networks than in interdependent Erdos Renyi networks. We 
are investigating whether the same effect appears in other opinion models of interdependent 
networks. The results will be presented in a future paper. 
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(a) t=2 




FIG. 1: Dynamics of the NCO model showing the approach to a stable state on a network with 
N = 9 nodes, (a) At t=0, five nodes are randomly assigned to be <r+ (empty circle), and the 
remaining four nodes are assigned with cr_ (solid circle). In the set comprising of node A and its 
4 neighbors (dashed box), node A is in a local minority opinion, while the remaining nodes are 
not. Thus at the end of this simulation step, node A is converted into <r_ opinion, (b) At t=l, in 
the set of nodes comprising node B and its 6 neighbors (dashed box), node B becomes in a local 
minority opinion, while the remaining nodes are not. Thus, node B is converted into <r_ at the 
end of simulation step t = 1. (c) At, t = 2, the nine nodes system reaches a stable state. 
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FIG. 2: Plots of s±, S2 and F of opinion cr + as a function of / with network si 
(a) ER networks with (k) = 4 and (b) SF networks with A = 2.5 and k } 
were done for 10 4 networks realizations. 
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FIG. 3: (a) Plots of n s as a function of s at criticality for ER networks with (k) = 4 and SF 
networks with A = 2.5 and k m [ n = 2. The dashed line is a guide to show that the slope obtained 
is t = 2.5. (b) Plots of Si as a function of N at criticality for ER networks with (k) = 4 and SF 
networks with A = 2.5 and k m [ Q = 2. The dashed lines are guides to show that the slope obtained 
is 6 ~ 2/3. All simulations were done for iV = 10000 and over 10 4 networks realizations. 
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FIG. 4: Plots of (a) s± and (b) S2 of opinion <r+ as a function of / for ER networks with different 
values of (k) for N = 10000. (c) s± and (d) S2 of opinion <r+ as a function of / for SF networks 
with different values of {k) for N = 10000 and A = 2.5. The solid lines are guides for the eyes. All 
simulations were done for 10 4 networks realizations. 
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FIG. 5: Plots of (a) si and (b) S2 of opinion o+ as a function of /, for ER networks with different 
values of W for (k) = 4. The solid lines are guides for the eyes. All simulations were done with 
N = 10000 and for 10 4 networks realizations. 
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FIG. 6: Schematic plot of the dynamics of the ICO model showing the approach to a stable state 
on a network with N = 9 nodes, (a) At t = 0, we have a stable state where opinion cr+ (open 
circle) and opinion cr_ (filled circle) coexist, (b) At t = 1, we change node 1 into a inflexible 
contrarian (filled square), which will hold <r_ opinion. Node 2 is now in a local minority opinion 
while the remaining nodes are not. Notice that node 1 is an inflexible contrarian and even if he is 
in the local minority he will not change his opinion. At the end of this simulation step, node 2 is 
converted into cr_ opinion, (c) At t = 2, node 3 is in a local minority opinion and therefore will 
be converted into cr_ opinion, (d) At t = 3, the system reaches a stable state where the system 
breaks into four disconnected clusters, one of them composed of six a_ nodes and the other three 
with one a + node. 
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FIG. 7: For ER networks with (k) = 4, plots of (a) s\ and (b) S2 as a function of / for different 
values of (ft under strategy I. Plots of (c) s± and (d) S2 as a function for / for different values of 
(ft under strategy II. The solid lines in S2 are guides for the eyes. All simulations were done with 
N = 10000 and 10 4 network realizations. 
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FIG. 8: F(k) as a function of k for different values of / for ER networks with (k) = 4 . We can see 
that as / increases F(k) — > 1 for smaller values of k. All simulations were done with N = 10000 
and 10 4 network realizations. 
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FIG. 9: For ER networks with (k) = 4, plots of (a) si (b)(top) S2 and (b) (bottom) AS\/ A<f) as a 
function of </> for different values of / under strategy I. Plot of (c) s\, (d)(top) S2 and (d) (bottom) 
A5i/A0 as a function of <fi for different values of / under strategy II. The solid lines in S2 and 
ASi/A<p are guides for the eyes. All simulations were done with TV" = 10000 and 10 4 network 
realizations. 
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FIG. 10: Plots of Tl g ELS db function of s for both strategies at f c {4>)- F° r the random strategy, 
fc(4>) = 0.36 and for the targeted strategy, f c {4>) = 0.45 when 4> = 0.2. The dashed line is a guide 
to show that the slope obtained is r ~ 5/2. All simulations were done with N = 10000 and 10 4 
network realizations. 
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(c) t=2 (d)t=3 (e)t=4 

^^^^^^^r ^^^^^l^r ^^^^^H^ 

FIG. 11: (Color online) Schematic plot of the dynamics of the NCO on coupled networks showing 
the approach to a stable state on a system of interdependent networks A and B with N = 6 nodes 
in each network. Open circles represent opinion <t+ and solid circles represent opinion <r_. The 
solid lines connecting nodes in each network are connectivity links within the networks, and the 
dashed lines connecting nodes from two networks are interdependent links. In the initial state, each 
node is randomly assigned with opinion er + or er_. (a) At t=0, opinion dynamics evolve within 
each single network. In networks A and B, nodes A\, A2, A3, A4 and B\ are in a local minority 
opinion within their network, so at the end of this time step these nodes will change their opinions. 
The remaining nodes will keep their opinions, (b) At t=l, the two networks interact through the 
interdependent links. Notice that the global majority opinion is <7+ now. Thus, the pairs A\ — B\ 
and A<z — B^ where two nodes have different opinions, will follow the global majority opinion. So 
at the end of this time step nodes A\ and A^ will change their opinions. The pair ^3 — .B3, remains 
as o"_ since both nodes share the same opinion, (c) At t=2, in network A, node ^3 is in a local 
minority, so at the end of the time step it will change opinion, (d) At t=3, the two networks 
interact through the interdependent links. Notice that the global majority opinion is still cr + , the 
only pair with different opinion is A3 — B3, s??at the end of this time step, -B3 will change its 
opinion, (d) At t=4, the interdependent networks reach a stable state and no more changes will 




FIG. 12: Plots of NCO on coupled ER networks, with (fc^) = (ks) = 4 and for each network 
N = 10000. (a) Plot of s\ of opinion <r + as a function of / for several values of q. As seen, when 
q = 0, the system undergoes a smooth second order phase transition (regular NCO model). As q 
increases until q = 0.5, it becomes a hybrid phase transition, which contains both a smooth second 
order type and a seemingly abrupt jump, i.e., s\ changes smoothly close to f c {q) and followed by 
a sharp jump at / = 0.5. When q is further increased (q > 0.5), the smooth phase transition 
disappears, the system undergoes a pure abrupt phase transition. In the inset of (a) we plot f c 
as a function of q. (b) Plot of S2 as a function of / for different values of q. As seen, when q 
increases, the peaks of the S2, which characterize the critical threshold value of the second order 
phase transition, shift to the right. We can see that beyond q = 0.5, the peak of S2 disappears, 
which indicates that there is no second order phase transition, (c) Plot of the change of S\ around 
/ = 0.5, AS*!, as a function of system size N for different values of q > 0. In the inset of (c), we 
plot AS\/N as a function of N for different values of q. The linear relationship between AS\ and 
N suggests that for q > 0, around / = 0.5, there exists a discontinuous transition, (d) Plot of the 
number of cascading steps (NOI) of the networks as a function of / for different values of q. As 
seen there are two peaks for each q value for q < 0.5, and for q > 0.5 there is only one peak at 
/ = 0.5. In the inset of (d), plot of the locatiorg4f NOI peak as a function of q. The solid lines are 
guides for the eyes. 
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FIG. 13: Plot of the number of cascading steps between the networks (ER with (Ua) = (^s) 
(k) = 4) at / = 0.5 as a function of system size for different values of q in log-log scale and in lo 
linear scale in the inset respectively. The dashed lines are the power law and logarithmic fittin 
respectively. 
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FIG. 14: Plot of n s as a function of s at criticality, f c , for different values of q for the NCO model 
on coupled ER networks, with (kji) = {ks) = {k) = 4 and N = 10000. For each network the results 
are averaged over 10 4 realizations. As q increases, n s losses the power law shape indicating that 
the second order phase transition is lost. The dashed line is a guide to show the slope r = 2.5. 
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FIG. 15: Study of NCO model on coupled SF networks, with k min = 2, A = 2.5 and N = 10000 for 
each network, (a) Plot of s± of opinion cr + as a function / for different values of q. (b) Plot of S2 
as a function of / for different values of q. (c) Plot of the number of cascading steps, NOI, of the 
coupled networks as a function of / for different values of q. The solid lines are guides for the eyes. 
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FIG. 16: Plot of n s as a function of s at criticality for different value of q for the NCO model 
on coupled SF networks, with fc m ; n = 2, A = 2.5 and N = 10000. For each case the results are 
averaged over 10 4 realizations. As q increases, n s losses the power law shape indicating that the 
second order phase transition is lost. The dashed line is a guide to show the slope r = 2.5. 



38 



